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ABSTRACT 

This paper presents the second release of Pico (Parameters for the Impatient COsmologist). Pico 
is a general purpose machine learning code which we have applied to computing the CMB power 
spectra and the WMAP likelihood. For this release, we have made improvements to the algorithm as 
well as the data sets used to train Pico, leading to a significant improvement in accuracy. For the 9 
parameter nonflat case presented here Pico can on average compute the TT, TE and EE spectra to 
better than 1% of cosmic standard deviation for nearly all £ values over a large region of parameter 
space. Performing a cosmological parameter analysis of current CMB and large scale structure data, we 
show that these power spectra give very accurate 1 and 2 dimensional parameter posteriors. We have 
extended Pico to allow computation of the tensor power spectrum and the matter transfer function. 
Pico runs about 1500 times faster than CAMB at the default accuracy and about 250, 000 times faster 
at high accuracy. Training Pico can be done using massively parallel computing resources, including 
distributed computing projects such as Cosmology@Home. On the homepage for Pico, located at 
http://cosmos.astro.uiuc.edu/pico, we provide new sets of regression coefficients and make the 
training code available for public use. 

Subject headings: cosmic microwave background — cosmology: observations — methods: numerical 



1. INTRODUCTION 

Given the quantity of data available from current ex- 
periments such as WMAP and SDSS as well as the 
prospects for the next generation Planck and DES ex- 
periments on the horizon, there is growing need for cos- 
mologists to develop tools that can accurately inter- 
pret the flood of data these experiments will gather. A 
key component in all such analysis is the exploration 
of the posterior density of the cosmological parameters 
given the available data. This allows us to constrain 
and test theoretical models of the Universe. A ma- 
jor computational hurdle in this procedure is the abil- 
ity to quickly and accurately compute the power spec- 
trum of CMB fluctuations and the matter transfer func- 
tion. This is accomplished u sing codes such as CMBFast 
(ISeliak fc Zaldarriagal Il996l ) and CAMB (jLewis et all 
l2000r ) that evolve the Boltzmann equation for the vari- 
ous constituents of the Universe. Using the default accu- 
racy settings in CAMB the calculation of a single power 
spectrum takes on the order of a minute. At higher 
settings, as may be required by upcoming CMB experi- 
ments, the computational time can jump to several hours 
on a modern desktop. Furthermore, computation of con- 
straints based on the data requires evaluating the power 
spectrum at O (lO^ — 10^) models. Decreasing the time 
required to calculate the CMB power spectrum, while 
maintaining sub-cosmic variance accuracy, will play an 
important part of turning raw data into quantitative in- 
formation about the history and structure of the Uni- 
verse. 

Previous codes aimed at speeding up power spec- 
trum computations such as DASH (jKaplinghat et all 
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I2002D and CMBwarp (jJimenez et al.l I2004D have at- 
tempted to reproduce the computation of CMBFast or 
CAMB. More recently, motivated to develop a code 
that was both faster than DASH and more accurate 
than CMB warp, we released a ma chine learning code 
called Pico (|Fendt fc WandeltJ[200l . It uses a training 
set of power spectra from CAMB to fit several multi- 
variate polynomials as a function of the input parame- 
ters. Along with accurately fitting the power spectra, 
we also found that Pico was able to directly fit the 
WM AP likelihood ^Spergel et al.l 120071 : iHinshaw etall 
120071: iPage et al.1l2007f). Thi s was done previously by 
CMBFit (jSandvik et al1l2004( l which used a similar idea 
of fitting the likelihood with a polynomial in the cosmo- 
logical parameters. By replacing CAMB and the WMAP 
likelihood code with Pico we demonstrated that it can 
quickly explore the parameter space and give nearly iden- 
tical posteriors. Since the first release of Pico, Auld et 
al. have applied a neural network code called CosmoNet 
to flat (Auld et al. 2006) and nonflat (|Auld et al.ll2007| ) 
models. Also Habib et al. have demonstra ted that a 
Gaus sian process model introduced in (Heitmann et al] 
[200l can predict the CMB power spe ctrum for flat mod - 
els based on a very small training set (jHabib et al.ll2007f l . 
Here we discuss some improvements to increase the ac- 
curacy of Pico and extend its to application to more gen- 
eral cosmological models. Further, we describe a training 
method that avoids serial runs of CAMB and is designed 
to leverage access to massively parallel and even dis- 
tributed computing resources. For example, we demon- 
strate that Pico can be trained using the thousands of 
geographically distinct hosts that contribute to Cosmol- 
ogy@Home.^ 

This paper is organized as follows. Section [2] summa- 
rizes some of the improvements to the Pico algorithm 
and how training sets are generated. In section [3] we 

http://www.cosmologyathome.org 
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demonstrate the performance of Pico in computing the 
power spectrum, matter transfer function and WMAP 
Hkehhood for a 9 parameter nonflat model. We show 
that Pico can be used to quickly explore the parameter 
posterior for this model while accurately reproducing the 
1 and 2 dimensional marginalized distributions. Lastly 
we summarize and conclude in section 

2. THE ALGORITHM 

2.1. Overview 

Given a training set of Cosmological parameters x and 
CMB power spectra y, Pico models the function y = 
/ (x) in 3 steps. First it compresses the power spectra 
using a Karhunen-Loev e (|Karhunenl Il946t iLoevd flOSSl : 
iTegmark fc Bunnlll995( l technique. For Cax = 3000, the 
temperature and polarization spectra due to scalar and 
tensor perturbations can be compressed to ^ 180 total 
numbers. Next the algorithm clusters the input param- 
eters into non-overlap ping regions. This is done us ing 
a fc-means algorithm (jMacQueenl 119671 : iKirbvl 1200 iD or 
by choosing hyperplanes to manually partition the space. 
Lastly, Pico models the function by fitting a least-squares 
polynomial to the compressed power spectra over each 
cluster. De tails of the algori t hm c an be found in the 
appendix of iFendt fc Wandeltl ()2007f ). 

2.2. Improved Fitting 

The first change to the algorithm is that it now uses the 
LAPACK libraries ^ to perform matrix decomposition. 
LAPACK makes the algorithm more stable, allowing the 
use of higher order polynomials. This gives Pico better 
fits at low e (< 200) and high £ (> 1500). This also makes 
clustering significantly less important but still useful for 
improving the low £ accuracy as well as improving the 
fit to the WMAP likelihood directly. In particular, we 
have found that partitioning the training set along values 
of constant ilk improves the low £ accuracy in the tem- 
perature power spectrum and matter transfer function, 
while partitioning in r, the reionization depth, gives a 
large improvement in computing the polarization power 
spectrum at low £. 

2.3. Decreasing Numerical Noise 

For nonflat models, fitting algorithms are hindered by 
the fact that at the default accuracy the power spec- 
tra computed by CAMB are numerically noisy. This is 
demonstrated in Figure [Tl where we have plotted the 
power spectrum at the default and high accuracy lev- 
els for various ^-values as a function of a single param- 
eter (r2b^^) for a nonflat model. Since the power spec- 
trum is not a numerically smooth function of the cos- 
mological parameters, Pico is limited in its ability to fit 
the true spectra. Also note that the higher accuracy 
power spectrum does not always smooth over the lower 
accuracy case. To reduce this noise and limit the ef- 
fect of interpolation errors we have generated training 
sets with CAMB using the accuracy parameters set as: 
Accuracy=3, lAccuracy=3 and lS£mip=3. In Figure [2]we 
have plotted the error between the default and high ac- 
curacy power spectra from CAMB for 25 models around 
the peak of the WMAP likelihood. The left and right 

^ http://www.nctlib.org/lapack/ 
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Fig. 1. — The plot shows the value of the temperature spectrum 
as a function of Cl\^h^ at various ^-values for a nonflat cosmology. 
The red (-I-) points correspond to to the default CAMB accuracies 
(1,1,1) and the green (□) points correspond to higher accuracy 
settings (3,3,1). While at low £ the power spectrum is smooth, the 
default accuracy becomes numerically noisy at higher £. This is 
one reason adding Slj^ as a free parameter increases the difficulty 
in fitting the power spectrum. Also plotted, as a blue line, is 
the power spectrum computed by Pico trained on the 9 parameter 
nonflat model discussed in section [S] 
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Fig. 2. — The plots show the percent error between the TT (left) 
and EE (right) power spectrum computed by CAMB at the default 
accuracy compared to those computed at high accuracy . Also 
shown is the percent error between the power spectra computed 
by Pico and the high accuracy CAMB spectra. This test was done 
on 25 models all located within 25 log-likelihoods of the WMAP 
peak. 



plots show the percent error in the TT and EE spectra. 
Also shown is the error between the spectra computed by 
Pico and the high accuracy CAMB runs. Here Pico was 
trained on the 9 parameter model discussed in section [31 
The figures demonstrate that when trained on high accu- 
racy data Pico computes the power spectrum around the 
peak of the likelihood to better precision than CAMB at 
its default accuracy settings. 

Numerical noise in the power spectrum also leads to 
noise in the WMAP likelihood as shown in Figure [31 
Again, this increases the difficulty in fitting the likeli- 
hood with Pico. Just as the power spectrum, this is 
remedied by running CAMB at higher accuracy. While 
the level of noise introduced in the likelihood may not be 
of significant concern when analyzing current CMB data 
sets, it may represent an important hurdle to overcome 
for the next generation of experiments. Also evident from 
Figure [31 is that Pico provides a smooth approximation 
to the noisy likelihood. This is an important property 
for algorithms that require differe ntiating the l ikelihood 
such as Hamiltonian Monte Carlo ()Haiianll2007i) . 

2.4. Generating the Training Set 
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Accuracy = 3 □ 
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Fig. 3. — Value of — InLwMAP ^ function of Q^^h'^ for a non- 
flat cosmology. Note that this is near the peak of the likelihood 
in the full space. The red (+) points correspond to the default 
CAMB accuracies (1,1,1) and the green (□) points corresponding 
to using higher accuracy settings (3,3,1). The blue line is the value 
computed by Pico trained over the 9 parameter nonflat models dis- 
cussed in section [3] Note that using the default accuracy in CAMB 
gives a numerically noisy function, which can lead to variations of 1 
or more log-likelihoods, but Pico gives a smooth function through 
the high accuracy values. 



As in the first release of Pico, we generate the train- 
ing set of power spectra and matter transfer function by 
samphng uniformly from a large box. After training Pico 
to compute the power spectra, evaluation of the WMAP 
likelihood, which requires a few seconds, becomes the 
new bottleneck in parameter estimation. Another signif- 
icant speed up can be obtained by using Pico to directly 
fit the likelihood function. However, for the 9 parame- 
ter case we examine in the next section this is a difficult 
problem. Training Pico over a box in parameter space in- 
cludes regions that are many thousands of log- likelihoods 
from the peak and gives a very sparse sampling of the 
high likelihood region. In practice we are not interested 
in these areas of parameter space. Instead we aim to 
compute the likelihood very accurately around the peak 
of the distribution. This requires generating a training 
set for Pico that includes only the high likelihood region. 

A natural method of accomplishing this is to use the 
Metropolis Hastings algorithm to find points in the high 
likelihood region. This can be done efficiently by running 
CosmoMCLewis fc Bridle (20Q2,) and using Pico to com- 
pute the power spectra. To ensure we cover a sufficient 
volume the chains are run using only the WMAP data 
and at a higher temperature, meaning the log-likelihood 
is scaled by a constant factor allowing the chains to ex- 
plore a larger volume. This step is dominated by the 
time it takes to run the WMAP likelihood code. Lastly 
we run the samples through CAMB and the WMAP code 
to get the true likelihood which will be used to retrain 
Pico. This step is also quick as it can easily be run in 
parallel. It is useful to note that this procedure never 
requires running CAMB in serial making it an ideal ap- 
plication for distributed computing projects such as Cos- 
mology@Home. The training set can be further refined 
by pruning out data at low likelihood. In the following 
Pico is trained to compute the likelihood on points within 
25 log-fikelihoods of the WMAP peak. 



2.5. Polynomial Hierarchy 

The process of generating the training set outlined in 
section 12.41 has the added benefit of giving us a set of 
power spectra constrained around the peak of the like- 
lihood. We would like to make use of these points by 
adding them to the power spectra training set. However 
adding a large weighting of points to a small region of the 
box will have a negative effet on the accuracy of the algo- 
rithm outside this region. Instead we have implemented 
the ability to use a hierarchy of polynomials with Pico by 
separately training over the uniformly sampled points in 
the full box and over only the points in the constrained 
region. If Pico is given a set of input cosmological pa- 
rameters within this region it computes the power spec- 
tra based on a polynomial fit to this constrained region. 
For points outside this region, Pico defaults to using the 
polynomials fit over the full box. While we have found 
that using only a single set of polynomials trained on the 
box is sufficient for analysis of current experimental data, 
this will be a useful feature in the future when data from 
higher resolution experiments become available. 

3. RESULTS 

Here we demonstrate the performance of Pico for non- 
flat cosmologies with the dark energy equation of state, 
wde allowed to vary (but still constant for a given 
model). In this space Pico fits the power spectrum and 
likelihood as a function of 
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The following sections study the accuracy of Pico in 
computing the power spectra, matter transfer function, 
WMAP likelihood as well as its application to parameter 
estimation based on this 9 parameter model. 

3.1. Power Spectra and Matter Transfer Function 

In order to demonstrate Pico's accuracy and robust- 
ness we will test the algorithm for two cases. The first 
case implements the hierarchy method discussed in sec- 
tion [231 Here the training set is divided into two pieces. 
The first contains ~ 18000 samples generated uniformly 
from the box defined in Table [l] and the second set con- 
sists of the 15000 points constrained to 25 log- likelihoods 
from the peak of the WMAP likelihood. For this case the 
test set consists of ^ 2000 points taken from the latter 
training set. These points were removed from the train- 
ing set and not used to train Pico. 

The models in the training set were run through 
CAMB at accuracy settings (3,3,3) to compute the true 
power spectra and transfer functions. As the £ and k 
sampling used by CAMB is model dependent it is nec- 
essary to spline the power spectra and transfer function 
so that each is computed at the same £ or k value. The 

values were chosen to be those used by CAMB for flat 
models with lSainp= 3. For the transfer function we 
used a unform sampling in Infc. Pico was trained using 
gth order polynomials, requiring less than 30 minutes on 
a 2.4GHz desktop. 

Pico's performance on this test set is shown in Figure 
H The top 2 rows show the TT, TE and EE power spec- 
tra with the second row focusing on low £. Results for 
the BB spectra and matter transfer function are shown 
in the third row. The two lines in each plot represent the 
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mean error and the error bar that bounds 99% of the test 
set. For the power spectra the error is plotted in units of 
the cosmic standard deviation and for the matter trans- 
fer function the y-axis shows percent error. From the 
figure we see that over the volume of parameter space 
important for CMB parameter estimation Pico can com- 
pute the power spectra for 99% of models in the training 
set to better than 4% of cosmic standard deviation, with 
the mean error around 0.5%, over most ^-values. Even 
the worst fits to the EE power spectra, which occur just 
after the reionization bump, are only about 25% of the 
cosmic standard deviation. We also note that many of 
the models that are fit poorly have very low power in 
this region so no experiment should be sensitive to these 
errors. For the transfer function the 99% error bar is 
around 0.25%, and the mean at 0.02%, except at very 
low k. This should be sufficient for analysis of data from 
the next generation of large scale structure experiments. 
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TABLE 1 

Parameter bounds defining the box the training set was 

sampled from for the example in section \3\ this 
encompasses a volume of at least 3(t in each parameter 
around the wmap maximum likelihood. note that we 
also impose the prior that the corresponding hubble 
constant for each parameter point lie in the interval 
[30, 100] which excludes some regions inside the box. 

For the second test case the hierarchy method is not 
used and Pico is only trained on a uniform sample of 
points from the box in Table [T] For this case the test set 
consists of a uniform sample of ^ 2000 points from the 
same box. Pico's performance on this test set is shown in 
Figure [5l We include this case only to allow comparison 
with other codes. When Pico is used to explore the pa- 
rameter posterior based on CMB constraints, which is its 
main purpose, chains will rarely propose points outside 
the constrained volume used in the hierarchy method. 
Therefore Figure S] provides a better indicator of the 
types or error incurred by using Pico to compute the 
power spectra. The regression files on the Pico website 
implement the hierarchy method. 

3.2. WMAP Likelihood 

Next we test the computation of the WMAP likeli- 
hood from Pico for two cases. The first uses Pico to 
compute the power spectrum and then the WMAP code 
to compute the likelihood and the second uses Pico to 
directly compute the likelihood. The training set for the 
likelihood computation consists of ^ 15000 points gener- 
ated using the method described in section [^T^ Another 
^ 2000 points, generated using the same method, were 
used as a test set. The absolute error between the like- 
lihood is shown in Figure [51 The plot on the left shows 
the results of using Pico to compute the power spectrum 
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Fig. 4. — The above plots compare the performance of Pico with 
CAMB at high accuracy settings for 9 parameter nonflat models 
with uiQE 7^ 1. Pico was trained using the hierarchy method de- 
scribed in section [2]5] and the test set consists of 2000 points within 
25 log-likelihoods of the WMAP peak. The top two rows show the 
error compared with CAMB in units of cosmic standard deviation 
for the TT, TE and EE power spectra at high £ (top) and low £ 
(center). The bottom row shows the error in the BE spectra in 
units of the cosmic standard deviation and the percent error in the 
matter transfer function. The two lines on each plot denote the 
mean error and the error bar that bounds 99% of the test set. We 
note that much of the error at low £ in the EE spectra is due to 
the 1% of models with extremely low power over this range. These 
spectra are too small to detect even with Planck. 




Fig. 5. — The plots are same as those in Figure |4] except here 
Pico was trained and tested over models sampled uniformly from 
the box defined by Table [T] Even over this larger region Pico can 
compute the power spectrum in 99% of the test cases to better than 
5% of cosmic standard deviation over most £ and is never worse 
than 0.7 cosmic standard deviation. 

while the plot on the right shows the results of directly 
computing the likelihood. For the case of directly eval- 
uating the likehhood, Pico can compute about 90% of 
the test set better than 0.25 log-likelihoods. When only 
using Pico to compute the power spectrum the results 
are within 0.25 log-likelihoods for better than 99.5% of 
the models. The training set and test set were computed 
using high accuracy CAMB runs and version v2p2p2 of 
the WMAP hkelihood code.^ 

3.3. Parameter Estimation 
^ http://lambda.gsfc.nasa.gov 
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Fig. 6. — The plots show the absolute error when evaluating 
the WMAP likelihood with Pico. In the left plot Pico was used 
to compute the power spectrum which were fed into the WMAP 
code. The right plot shows the absolute error when Pico directly 
evaluates the likelihood. In both cases the likelihood is compared 
to the value of the WMAP code using high accuracy CAMB power 
spectra. 



To test the application of Pico to this 9 param- 
eter model we ran Markov chain s usin g CosmoMC 
with the W MAPdHinshaw et all 120071: iPage et al 



ACBAR dKuo et al.ri2004f) CBIdReadhead et al 



SDSS 



BoomerangfMontrov ct al.' 2006; Jones ct aL 

fcentini ct al. 2006), 2df(Cole et al. 2005), 
iclman- McCarthv et all l2006f) . and SNLS 
(|Astier et al., ,2005i) data sets. Chains were run for 3 
cases. The first uses CAMB and the WMAP likeli- 
hood code (CAMB-I-WMAP), the second uses Pico to 
compute the power spectra and transfer function but 
still uses the official likelihood codes (PICO-t-WMAP) 
and the third case uses Pico to compute the power 
spectra, transfer function and the WMAP likelihood 
(PICO). In third case we did not use Pico to fit the 2df, 
SDSS or SNLS likelihood codes. The 1-dimensional, 
marginalized posteriors for each of the 9 parameters 
are shown along the diagonal in Figure [T] The plots 
in the lower (upper) triangle in the figure compare the 
2 dimensional posteriors between the PICO-fWMAP 
(PICO) case and the CAMB+WMAP case. In all of 
the plots the CAMB-I-WMAP results are shown in 
red, the PICO+WMAP results in green and the PICO 
results in red. The lines in the 2D plots denote the 
68% and 99% contours. Using 6 chains run in parallel, 
the PICO-I-WMAP chains ran about 60 times faster 
than without Pico, requiring about 4 hours of wall clock 
time. Using Pico to compute the WMAP likelihood 
gave another factor of 2.5 decrease in CPU time giving 
a total speed up of ~ 150 over chains run without Pico. 
In all cases the chains finished with a Gelman-Rubin 
statistic less than 1.01. 

4. CONCLUSION 

This paper describes a major new release of Pico, a 
fast and accurate code for computing the CMB power 
spectrum, matter transfer function and the WMAP like- 
lihood. We noted the presence of numerical noise in 
CAMB at standard accuracy for nonflat models and its 
effect on the WMAP likelihood. To solve this problem 
we have generated training sets running CAMB at high 
accuracy settings. Also we have presented a method of 
generating a training set that finds the high likelihood 
region of parameter space without ever running CAMB 
in serial. This is especially useful for training Pico to fit 
the WMAP likelihood in large dimensional spaces. Fur- 
thermore, Pico can be trained separately on the power 
spectra in this smaller region of parameter space allowing 
even more accurate results around the peak of the like- 



lihood while still maintaining the ability to compute the 
power spectra over a large box in parameter space. The 
combination of these improvements, along with modifica- 
tions to the Pico algorithm, have increased its accuracy 
in computing the power spectrum and likelihood. Also 
we have extended Pico to compute the power spectrum 
due to tensor perturbations as well as the matter power 
spectrum. On the Pico homepage,^ we provide the new 
version of Pico and new sets of regression coefficients. 
We have also released the training code for Pico, allow- 
ing users to apply the algorithm to new classes of models 
and parameter sets. We expect that the accuracy and 
speed achieved by Pico will be useful for current and 
future CMB and large scale structure observations. Fur- 
thermore, we hope that the concept, embodied by Pico, 
of exploiting massively parallel computing resources to 
solve inherently serial numerical problems will find ap- 
plications beyond the immediate domain of cosmological 
parameter estimation. 

This work was partially funded by NSF grants 
AST 05-07676 and AST 07-08849, by NASA contract 
JPL1236748, by the National Computational Science Al- 
fiance under AST300029N, by the University of Illinois, 
by the Computational Science and Engineering Depart- 
ment at the University of Illinois and by a Friedrich Wil- 
helm Bessel research prize from the Alexan der von Hum- 
boldt foundation. We utilized the Teragrid llCatlett et al.l 
l2007f) Itanium 2 clusters at NCSA and at Argonne Na- 
tional Laboratory, as well as the Turing cluster in the 
Computational Science and Engineering Department at 
the University of Illinois at Urbana-Champaign. 

We thank the Max Planck Institute for Astrophysics 
for its hospitality while part of this work done. We also 
thank the users of the Cosmology@IIome^ project whose 
donated CPU hours helped make this work possible.^ 
In particular we would like to thank Scott Kruger, the 
administrator of Cosmology@IIome, as well as the users 
laurenu2, PoorBoy, [B^S]ralfi65, Mitchell, and Mike The 
Great as representatives of all Cosmology@Home partic- 
ipants. Lastly we would like to thank Nikita Sorokin for 
his work in designing the Pico homepage. 

Funding for the Sloan Digital Sky Survey (SDSS) 
has been provided by the Alfred P. Sloan Foundation, 
the Participating Institutions, the National Aeronautics 
and Space Administration, the National Science Foun- 
dation, the U.S. Department of Energy, the Japanese 
Monbukagakusho, and the Max Planck Society. The 
SDSS Web site is http://www.sdss.org/. The SDSS 
is managed by the Astrophysical Research Consortium 
(ARC) for the Participating Institutions. The Participat- 
ing Institutions are The University of Chicago, Fermilab, 
the Institute for Advanced Study, the Japan Participa- 
tion Group, The Johns Hopkins University, the Korean 
Scientist Group, Los Alamos National Laboratory, the 
Max-Planck-Institute for Astronomy (MPIA), the Max- 
Planck-Institute for Astrophysics (MPA), New Mexico 
State University, University of Pittsburgh, University 
of Portsmouth, Princeton University, the United States 
Naval Observatory, and the University of Washington. 



http 

8 http 

9 http 



/ /cosmos. astro. uiuc.edu/pico 
/ /www. cosmologyathome.org 
/ / www.cosmologyathome.org/top_users.php 



6 



0.022 0.026 0.085 0.11 1.037 1.052 0.048 0.161 -0.034 0.009 -1.26 -0.54 0.953 1.06 2.891 3.152 0.175 




0.022 0.026 0.085 0.11 1.037 1.052 0.048 0.161 -0.034 0.009 -1.26 -0.54 0.953 1.06 2.891 3.152 0.175 0.7 

Fig. 7. — The cosmological parameter posteriors using CAMB and Pico for 9 parameter nonflat models with uide ~1 based on the 
WMAP, ACBAR, CBI, Boomerang, SDSS, 2df and SNLS data sets. The red Unes are the result of using CAMB and the WMAP likelihood 
code. The green lines use Pico to compute the power spectrum but still uses the WMAP likelihood. The blue lines are the result from 
using Pico to compute the power spectrum and WMAP likelihood. The plots in the lower triangle show the 68% and 99% contours for the 
chains run using CAMB and PICO with the WMAP likelihood code. The upper triangle compares the CAMB chains to those using Pico 
to compute the power spectra and WMAP likelihood. 
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